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The D'yakonov-Perel' mechanism of spin relaxation is connected with the spin splitting of the 
electron dispersion curve in crystals lacking a center of symmetry. In a two-dimensional noncen- 
trosymmetric system, e.g. quantum well or heterojunction, the spin splitting is a linear function 
of k, at least for small values of k. We demonstrate that the spin relaxation time r s due to the 
spin splitting is controlled not only by momentum relaxation processes as widely accepted but also 
by electron-electron collisions which make no effect on the electron mobility. In order to calculate 
the time r s taking into account the electron-electron scattering we have solved the two-dimensional 
kinetic equation for the electron spin density matrix. We show how the theory can be extended to 
allow for degenerate distribution of the spin-polarized two-dimensional electron gas. 

I. INTRODUCTION 

In recent years the electron spin relaxation processes have received much attention in connection with various spin- 
tronics applications. For free elec trons in quantum well structures the following four mechanisms of spin decoherence 
are most important (see Refs. llld and references therein): the Elliot- Yafet, D'yakonov-Perel', Bir-Aronov-Pikus mech- 
anisms and electron spin-flip scattering by paramagnetic centers. This paper is devoted to the D'yakonov-Perel' spin 
relaxation (DPSR) in which case the spin relaxation time, r s , is related to the spin splitting and given by Ref. 

t- 1 oc (nl) T . (i) 

Here f2fc is the effective Larmor frequency describing the splitting of the electron spin dispersion branches, the angle 
brackets mean averaging over the electron energy distribution and r is a microscopic electron relaxation time. In a 
two-dimensional (2D) system laking a center of a symmetry, e.g. quantum well or heterojunction, the frequency fife 
is linear in k. In this case the time r was usually identified with the momentum relaxation time that determines the 
electron mobility Refs. 00000- In Ref. we have shown that the inverse time r _1 is determined not only by the 
momentum scattering rate but contains also an independent contribution from electron-electron collisions which make 
no effect on the electron mobility. Really, electron-electron collisions change the direction of k and ilfc and, therefore, 
they control the Dyakonov-Perel' spin relaxation exactly in the same way as any other scattering processes do. 

The paper is organized as follows. Section 2 contains the discussion of mechanisms of the spin splitting of the electron 
subbands in quantum wells and electron spin decoherence due to spin splitting, Sec. 3 describes the D'yakonov-Perel' 
spin relaxation mechanism in terms of the kinetic theory, in Sec. 4 we analyze the electron-electron collision integral 
for the spin-polarized electron gas, and Sec. 5 describes the solution of the kinetic equation with the account of 
electron-electron collisions. 

II. SPIN SPLITTING OF ELECTRON SUBBANDS 

In the parabolic approximation the effective electron Hamiltonian in the zinc-blende-based (OOl)-grown quantum 
well (QW) can be written as 

h 2 k 2 

U = + Piidyky - a x k x ) + f3 2 (cr x k y - a y k x ) , (2) 

where as|| [100] , y||[010], /?i,2 are constants, m is the in-plane effective mass, ai are the Pauli matrices and k n are 
the components of the 2D electron wavevector, k. The terms with higher powers in k are not considered here. In 
the symmetrical QWs (the D 2 d point group) the linear-fc spin-dependent term proportional to (3i is allowed only, 
it is called the bulk-inversion asymmetry (BIA) term. In asymmetrical QW structures characterized by the C 2v 
point-group symmetry, there exists another spin-dependent contribution proportional to /?2 and called the structure- 
inversion asymmetry (SIA) term or the Rashba term Refs. ITollTll (see also Ref. ^3 an d references therein). The 
structure asymmetry can be related with non-equivalent normal and inverted interfaces, external or built-in electric 
fields, compositionally stepped QWs etc. 

It is convenient to introduce the Cartesian coordinates x'||[110], j/'||[110], z||[001] which allow to write a sum of the 
BIA and SIA terms in the form 

Wci(fc) = (P-Vx'ky + l3+Oy,k xl ) /2 , (3) 
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where f3± = 2(/?2 ± Pi)- The effective Larmor frequency fl k is defined by Tt c i — (h/2)ttk ■ cr, its components are 

Qk,x> = P-k y </h , Cl k , y > = (3 + k x >/h , Cl h ,z = 0. (4) 

The spin splitting at k — (k x i,k y i) is hflk = \jp\k\i + P-k 2 ,. If only one of the linear-fe terms, BIA or SIA, is 

nonzero then P 2 _ — P\ = p 2 and the splitting Ml k = pk is angular independent. 

Consider an electron gas occupying the lowest conduction subband el and assume that, at the moment t = 0, the 
electrons are spin-polarized in the same direction along, say, the growth axis z. Owing to the spin-splitting of the 
electron subband the electron spin in the state with the wavevector k precesses around the axis fife which lies in the 
interface plane, see Eq. @. In case of the large splitting, \flk\ S> 1/r, the spin of electron in the state k will oscillate 
with time as s z (t) = s z (0) cosfifci. In the case of being isotropic in the QW plane, the spin polarized electrons 
which occupy the circle of the fixed radius in the k space show the similar oscillatory behavior for their s z component. 
However, if P_ ^ 0i_ and/or the electrons occupy states with different values of |fc|, the scatter in results in a fast 
non-exponential decay of s z . 

D'yakonov and Perel' Ref. were the first to show that the processes of electron-momentum scattering result in 
the slowing off the spin decohcrcncc caused by the spin splitting. In the collision-dominated ( "motional- narrowing" ) 
limit | fife | <C 1/t, this results in an exponential decay of the average spin: s z (t) = s z (0) exp (— t/r a ), where r s is g iven 
by Eq. (QJ. The dimensionless coefficient in Eq. can be obtained from the solution of kinetic equation Refs. I3l4l 
This approach is valid as long as the energy relaxation time r e is short as compared to t s . 



III. DPSR IN TERMS OF THE KINETIC THEORY 



In the frame of kinetic theory, the electron distribution in the wavevector and spin spaces is described by a 2 x 2 
spin-density matrix 



Pk = fk + Sk ■ c- 



(5) 



Here f k = Tr(/5fc/2) is the average occupation of the two spin states with the wavevector fc, or distribution function 
of electrons in the fc-space, and the average spin in the k state is s k = Tr[/5fc(cr/2)]. If we neglect the spin splitting 
then, for arbitrary degeneracy of an electron gas with non-equilibrium spin-state occupation but equilibrium energy 
distribution within each spin branch, the electron spin-density matrix can be presented as 



P°k 



exp 



Ek — M ~ P (vos 



k B T 



(6) 



where E k = h 2 k 2 /2m, ks is the Boltzmann constant, T is the temperature, o s is the unit vector in the spin polarization 
direction, fi± = p, ± jl are the effective Fermi energies for electrons with the spin component 1/2 or —1/2 along o s so 
that the energy distribution functions of electrons with the spin ±1/2 are given respectively by 



fk,± = 



exp 



Ek — M± 
k B T 
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Note that Eq. © can be rewritten in the equivalent form Ref. |3| 

1 



P°k 



[/*,+ + fk - + (fk,+ - fk-) (o-o s 



The densities n± of 2D electrons with a particular spin can be related with the effective Fermi energies by 



n± 



2-irh' 



k B T\n(l 



H±/k B T 



(7) 



If the spin splitting is non-zero but small compared to fi/r, the distribution function Tr[pfc/2] = does not 
change, whereas the spin vector obtains a correction 6sk = Sfc — proportional to the spin splitting. Therefore, the 
spin-density matrix may be presented as 



Pk = P k 



5s 



k ■ cr 



(8) 
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The quantum kinetic equation for the spin-density matrix taking into account the electron-electron collisions has 
the form 

^ + ^{H cl (k),p k } + Q k {p} = 1 (9) 

where [P,R] — PR — RP, H c i(k) is the linear-fc Hamiltonian and the third term in the left-hand side is the 
collision integral or the scattering rate, in this equation it is a 2x2 matrix. It follows from Eq. @ that the pseudovector 
kinetic equation for s k can be written as 

^t + n k x s k + Q k {s} = 0, (10) 

where Q k {s, /} = (l/2)Tr(crQ k {p}). In the following we ignore spin flips under scattering. Then, say, for the elastic 
scattering one has 

Q k {F} = Y,W klk {F k ~F k ,) 1 

k' 

where W k ' k is the probability rate for the electron transition from the state fe to fc'. The collision integral for electron- 
electron scattering is considered in the next section. Here it suffices to note that, for the distribution p k , the collision 
integral vanishes identically. This integral also vanishes after the summation over fc which allows, in particular, to 
derive from Eq. (|10|l the following equation of balance for the total average spin So = o s (n + — n_): 

^+Y / n k x6s k ^0. (11) 

k 

The angular dependence of the non-equilibrium correction 5 s k is the linear combination of cos = k x i / k and 
sin<I>fc = kyi/k, where <S> k is the angle between fc and the axis x' . Retaining in the kinetic equation IjlUfl terms 
proportional to the first angular harmonics we obtain the equation for Ss k with the inhomogeneous term linear in 
So. Then one can substitute the solution in the second term of Eq. (TTJ. The final result is that the tensor of inverse 
spin relaxation times, is diagonal in the coordinate system x',y', z and given by Refs. 

where t p is the momentum relaxation time. If among the two contributions, BIA and SIA, to the spin splitting one 
is dominant and = |/?-|, the spin relaxation times are interconnected by Ref. U 



= 2t. 



zz 



Interplay between the BIA and SIA contributions can lead to a giant spin relaxation anisotropy Ref.ll3l In particular, 
if these contributions coincide, Pi = /3 2 , so that /3_ = one has t*, x , = t| z and t^, , = oo. In the case f3\ = —f3 2 the 
coefficient j3 + = 0, the time r*, x , is infinite and r*, ; coincides with r| z . 



IV. ELECTRON-ELECTRON COLLISIONS IN QWS 

Here we will write the electron-electron collision term Q k {p} in the general case of arbitrary spin-density matrix 
p k (in particular, arbitrary degeneracy and arbitrary distribution of spin in the fc-space). For this purpose we remind 
that the matrix element of the Coulomb scattering fc, s k + fc', s k i — > p, s p + p', s p < is given by 

M(p, s p ;p',s p >\k,s k ;k',s k ,) = V k - P S Sp<Sk 6 Sp , <Sh , - V k - P > S Sj)tSh ,S Sp , lSk , (13) 

where s k , s k i... — ±1/2, V q is a Fourier transform of the 2D Coulomb potential of the electron-electron interaction 

2ne 2 , , 

= aegE ' ^ 

e is the elementary charge, as is the dielectric constant, and E is the sample area in the interface plane; in the following 
we set E = 1. Equation i|13|) takes into account both the direct and exchange Coulomb interaction. 
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In order to present Q k {p} in a compact form we introduce the 2x2 unit matrix /W and Pauli matrices 
(a = x, y, z) for the spin coordinates s pi s k and a similar set of four matrices, 1^ and o~a\ for the spin coordinates 



s p > , s k , . One can check that Eq. allows the following matrix representation 

M = AI^J^+Ba^-a^ , (15) 

A = V h - P - X - V k - P , , B = — i V k -p' . (16) 
Now the collision term for the electron spin-density matrix can be presented as 

QfcM = ^r $k +k >, P+P > S(E k + E k , - E p - E p ,)Tr 2 G(p,p';k,k') , (17) 

h'pp' 

G(p,p';k,k') = (18) 



- MpP$M(I (1) - ftW 2) P k 2 )) Pi iy )(l {2) $)M^$M . 

Here the spin-density matrices p^\k) = I^'fk + ■ s k , p^(k') — + <J ^ ' s k' c ^ c -i the symbol Tr 2 means 

the trace over the spin variable 2. After the trace is found the index 1 in Ti2G(p,p'; k, k') can be omitted. In order 
to derive Eqs. (|T7Jl . (fT5|) we used the standard diagram technique. 

Instead of equation @ for the spin-density matrix one can use a scalar equation for the distribution function fk in 
the form 

^- + Qk{f,s} = (19) 

and a equation for the spin-distribution vectorfunction as 

ds l» 

-f + n k x s k + Q k {s,f} = 0, (20) 
at 

where 

Qk{/,s} = ^Tr 1 [Q fe {/5}] (21) 
= ^ J2 <W,p+p' 5{E k +E k , -E p - E p ,) Tr 1 Tr 2 [G(p,p';fc,fc')] . 

k'pp' 

Q k {sJ} = ±Tr 1 l<T^Q k {p}} (22) 
= ^ 5 k+k', P +p> 5{E k + E k , -E v - Ep,) Tr^a^ G(p,p'; k,k')] . 

k'pp' 

For the analysis of the general equations 0, (|1 T|) . 1|21|1 . (122I) we consider below few particular cases. 

Spin-unpolarized electrons. In this case s k = and the spin-density matrix reduces to a product of the unit 2x2 
matrix and the distribution function f k . Taking into account that 

M 2 = (A 2 + 3B 2 ) I^I (2) + 2B{A - B) cr (1) ■ a (2) 

we come to the conventional collision term 

Qk{/} = W. p+p' 5(E k + E k , -E v - Ep,) {A 2 + 3B 2 ) (23) 

k'pp' 
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x[/ fc / fc '(l - /p)(l - f P >) - /p/ P '(l - /fc)(l ~ /*)] ■ 
Note that according to Eq. i|16|) one has 

A 2 + 3B 2 = V 2 _ p + V 2 _ p , - V k _ p V k _ p , 

which is one-forth of the function R introduced in Eq. (2.4b) in Ref. and the above collision term agrees with the 
equation (2.4a) in the cited paper. It is worth to note that in the sum (|2lfl) V k _ p + V£_ p , can be replaced by 2V J ?_ p . 

Electrons polarized along the same axis. By using the coordinate system with z parallel to the electron spin 
polarization one has s kyX = s ktV = and the spin-density matrix is a diagonal matrix with the diagonal components 
f h>s (s = ±1/2). It follows then that the products (k)p^ (k') and [/W - p {1) {p)][I {2) ~ P {2] {p')] are diagonal as 
well. We can take into account the spin conservation 

M(p, s 3 ;p',Si\k,sx;k',S2) oc S S3+SiiSl+S2 

and use the identity 

M(p, s 3 ;p', s 4 \k, si; k' , s 2 )M(p, s 3 ;p', s 4 |fe, -sy, k' , s 2 ) = . 
This allows, in agreement with Ref. Il7l rewrite the collision term for f k >s as 

-j- ^ ^ 5k+k>, P + P ' 5(E k + E k > - E p - E p >)W(p,s 1 ;p',s 2 \k,s; k' , s') 

k'pp' s'siS2 

X[fk,sfk',s'0- - /p,si)(l - fp'.s 2 ) - /p,si/p',s 2 (l _ /fe,s)(l - /fe',s')] • 

Here 

W{p,s;p',s\k,s;k',s) = (A + B) 2 = (V k - P - Vk-p'f , 
W(p, s; p', -s\k, s; k' , -a) = (A - B) 2 = V 2 _ p , 
W{p,-s;p',s\k,s;k',-s) = (2B) 2 = V 2 _ p , , 

and other values of W with s\ + s 2 ^ s + s' are zero. 

Low electron polarization. If the average electron spin s k is small as compared with the occupation probability f k 
then, in the equation for f k , one can ignore the spin polarization at all and use Eq. (|23|l . while, in the equation for 
Sfc, one can retain in the collision term only the contribution linear in s k . The linearized collision term is given by 

Qk{sJ} = Y Sk+k',p+ P >S(E k + E k ,-E p -E p ,) (24) 

k' .p.p' 

x [{Vk-p + Vk-p- ~ V k - p V k - p/ )s k F(k';p,p') - V k - p V k - p ,s k ,F(k;p,p') 
-V k - p {V k - p - V k - p >)s p F(p';k,k') - V k - p ,{Vk- P > - V k - p )s p ,F(p;k,k')] , 

where 

F(k x ;k 2 , k 3 ) = f kl (1 - / fe2 )(l - / fe3 ) + (1 - f kl )f k J k3 = f kl (1 - f k2 - f k3 ) + f k J k3 . 
Equation 124|l can be transformed to 

Qk{s, f} = Y Yl 6 k+k',p+ P - 5{E k + E k , -E p - Ep,) (25) 

h' ,p,p' 

x {2^ fe 2 _ p [s k F(k';p,p') - s p F{p'-k,k')] - 



-Vk-pVk-p> [s k F(k l -p,p l ) + s k ,F{k l -p,p l )-2s p F{p l -k,k 1 )]}- 
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Here the term proportional to 2V k _ p is due to the direct Coulomb interaction whereas the term proportional to 
Vfe_pVfc_ p / comes from the exchange interaction. 

Non-degenerate 2D electron gas. In this case the function F(ki\ &2, k^) reduces to f kl and the collision terms are 
as follows 

2ir 

Qk{f, s} = yY, S k+k', P+P ' S(E k + E k , -E p - E„) (26) 

k'pp' 

X [( 2 ^fe- P _ Vk- p Vk-p')(fkfk' - fpfpi) - Vk-pVk-p' (sfc • Sfe/ - s p ■ s P ')] , 

Q k {s, /} = y E W'.p+p' 8(E k + E k , -E v - E p .) (27) 

k'pp' 

X [ 2 ^fc-p( s fc/fc' _ s p/p') _ Vk- P Vk-p'{skfk' + Sfc'/fc - 2s p /p')] . 

Neglecting the exchange interaction given by the term proportional to the V k - p V k - p > the scattering rate Q k {f, s} is 
independent on the electron spin distribution. Both Qk{f, s} and Q k {s,f} take a simple form 

Qk{f,s} = ^2 Wpp^kk'ifkfk' - fpfp') , Qk{sJ } = E Wp p , Mk >{s k f k , - Sp/p-) , (28) 

k'pp' k'pp' 

where W pp t k y is the probability rate for the scattering of a pair of electrons from the k, k' states to the p, p' states 

2ir 

W pp , ykk > = —S k+k > , p+p > S(E k + Ey - E p - E p /) 2 V k _ p , 
an additional factor of 2 takes into account the double degeneracy of the electronic states. 

V. ELECTRON-ELECTRON SCATTERING TIME CONTROLLING THE DPSR 

In what follows we consider only a non-degenerate 2D electron gas in which case the zero-approximation spin-density 
matrix (JSJ) can be written as 

pi = /°(1 + 2s ■ a) , (29) 

where is the Boltzmann distribution function, and s is the average spin per electron, S /(n + + n_). The nonequi- 
librium correction Ss k satisfies the equation 

n k x (2f°s )+Q k {6s,f°} = Q. (30) 

It should be noted that in the 2D case the collision term l|27[l does not allow the quasi-elastic and relaxation time 
approximations and Eq. (|3U[I must be solved directly. 

The a-component of the vector product £l k x Sq can be represented as (tt k x So) a = Aapyhpso^, where the third-rank 
tensor A in general case of both BIA and SIA linear-fc terms has four nonzero components 

A-xxz = —A-zxx — P+/Tl, A-yyz = — h-zyy = P-/fl ■ 

The function (l/kf3)Q k {kpF k , /°} is independent of the azimuthal angle (here F k is an arbitrary function of 
k = \k\) as the operator Q k {Ss a , f } conserves the angular distribution in the k space. In such case the solution may 
be written as follows 

6s a (k) = -A Q/37 ^s 07 k T r* ee e^ T v(K) . (31) 

Here we introduced the dimensionless wavevector K = fe/fc-r , fc-r = (2mfcsT/?i 2 ) 1 ' /2 , 

* frk B Tee 2 

= ~^N~ ' (32) 
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N = n + + n- and v(K) satisfies the equation 

Ke- R2 = J d 2 K' J d 2 PW PP ,, KK , (v{K)er K ' 2 - cos 9 v{P) e - p ' 2 ^ , (33) 

where is the angle between K and P, P' = K + K' P, 

H'pf.kk' = |K * 2 S(K 2 + K' 2 -P 2 - P 12 ) . 

Inserting Eq. I|31|l into Eq. I jllj l one obtains after summation over k the principal values of the tensor of reciprocal 
spin relaxation times 

1 W^M%, '_>_+>_. (34, 



^ / T y'y' V ^ / T zz T x'x' T y'y' 

The time r which controls the spin relaxation is given by 

oc 

t = r* ee I, I=\J v(K)K 2 dK . (35) 
o 

The parameter r* is present also in the ee-scattering time which determines the rate of energy exchange between 2D 
electrons Ref. Ha 

The function v(K) was expanded in series using a basis set l n (e) — \/2 exp (—e)L n (2e), where L n (e) are the Laguerre 
polynomials and e = K 2 . The expansion was substituted into the right-hand side of Eq. H33|) . the integration has 
been performed by Monte-Carlo method. The problem was reduced to a set of linear inhomogeneous equations for 
the expansion coefficients of v{K). The resulting value of I in Eq. H35[) was found to be w 0.027. Allowance for the 
exchange interaction leads to a slight increase of this value to 0.028. 

Comparison with ionized impurities scattering. If the ionized impurities of the same concentration N lie inside 
the 2D layer then the corresponding transport time is given by T tr = (2/7r 2 )r* e (see Ref. El). The spin relaxation 
time, r] 4 , controlled by scattering by the ionized impurities is given by Eq. (|34|l where r is changed by Tt r /2, thus 
the ratio of the spin relaxation time, r| e , governed by electron-electron collisions and the time rf is w 3.6, i.e. the 
clastic scattering by impurities is less efficient. If the doped layer is separated from the quantum well by a spacer, the 
influence of the Coulomb potential of ionized impurities on t s is reduced and eventually can be neglected. 

The above result was obtained for the 2D Coulomb potential U 2D {p) = e 2 /eep, where p is the distance between the 
electrons in the interface plane. It is this potential that leads to the Fourier transform given by Eq. I|14|). In order 
to analyze the role of the quasi-2D character of the electron wave function confined in a QW of the finite thickness a 
one can replace U 2D (p) by the effective potential obtained by averaging the three-dimensional Coulomb potential as 
follows 

e 2 f f 



where ip e i (z) is the electron envelope function at the lowest conduction subband el. The straightforward calculation 
shows Ref. 18 that for the conduction band offset V c S> h 2 /ma 2 and fc^a < 1 a value of the time r = r* e J increases 
with widening the QW as the electron-electron interaction becomes weaker, but an order of magnitude of I and r 
remains the same as in the exact 2D case. 



Conclusion. Future Work 



We have shown that electron-electron collisions control the D'yakonov-Perel' spin relaxation in the same way as any 
scattering processes do. The electron-electron collision integral has been derived for an arbitrary degeneracy and spin 
distribution of the 2D electron gas. The calculations has been performed for non-degenerate 2D and quasi-2D electrons 
confined in a quantum well. Calculations of the spin relaxation time controlled by electron-electron collisions in bulk 
semiconductors are in progress. An important next step is the extension of the calculations from non-degenerate to 
degenerate spin-polarized electron gas. 

In agreement with our theory, the latest optical spin-dynamic measurements in an n-doped 
GaAs/AlGaAs QW of high mobility give an experimental evidence of the electron-electron scattering to randomize 
the electron spin precession Ref. [l!J. 
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Note in conclusion that the time r = T* e I can be related physically with the momentum relaxation time of an 
electron by equilibrium holes of the density N if the electron and hole effective masses are assumed to coincide. In 
the electron collisions with holes the directed electron momentum is transferred to the hole gas and decays within the 
time ~ r. The work was supported by RFBR, and by the programs of the Russian Ministry of Sci. and the Presidium 
of the Russian Academy of Sci. 
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